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NUCLEAR CASCADES IN ELECTROMAGNETIC SHOWERS PRODUCED 
BY PRIMARY GAMMA-QUANTA IN THE ATMOSPHERE. 

Danilova T.V., Erlykin A.D. , Mironov A.V., Tukiah E.I. 

P.N. Lebedev Physical Institute, Leninski pr. 53 , 

Moscow, 117924. USSR. 

Distributions have been calculated for the number of 
electrons N e , number of muons with the energy above 5 GeV 
% and the energy of hadron component E k in electromagnetic 
showers, produced by primary gamma-quanta with energies 
10‘ - lO^eV, incident on the atmosphere at zenith angles 
0 ^ 30° and observed at the mountain level 700 g/cm 2 -. The 
mean number of nuclear interactions of photons with the 
energy above 5 GeV is about 0.3 per each TeV of the prima- 
ry energy and nuclear cascades take out in average about 
2% of the total shower energy. The mean number of 5 GeV 
muons for the electromagnetiqjshower is (2-5)% from the num- 
ber of muons in EAS with the same number of electrons at 
the observation level. Similar value for the total energy 
of hadron component is also (2-5)%. Nm. and N e values as well 
as E k and N a don't correlate at the fixed primary energy EZ. 
Between N„ and E k there is a positive correlation at the 
given E®. 

Simulation of nuclear cascades arising when the electro- 
magnetic shower propagates through the atmosphere has been 
undertaken in connection with the works on high energy gam- 
ma-astronomy carried out at the Tien-Shan EAS array /I/. 
Similar study is planned at ANI array /2/. The difference 
between the present set of calculations and that carried 
out in pioneer works of Wdowczyk /3/ is in their orienta- 
tion onto Tien-Shan and ANI arrays, i.e. the observation 
of showers at the mountain level 700 g/cm 2 , higher muon 
threshold energy - 5 GeV and the usage of hadron component 
as an additional criterium for the selection of pure elec- 
tromagnetic showers. 

Electromagnetic showers were simulated for primary gam- 
ma-quanta with energies 10^,10^,10^ and lO 1 ^ eV, incident 
on the atmosphere isotropically in the zenith angle inter- 
val 0 ^ 30®. The total number of electrons at the level 
2 

700 g/cm was determined after the sampling of the first 
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interaction point and of the value cos 0 by means of 
approximation formulae of Greisen /4/ or Kliraakov-Pavlov 
/ 5 /. Then the development of the electromagnetic shower 
in the atmosphere was followed by steps A z-10 g/cm • 

The probability for photons P (E® } > 5 GeV, z) to have pho- 
tonuclear interaction in ide the step was estimated accor- 
ding to formulae of the Approximation A of the cascade the- 
ory: W»1-exp(- P(B®, >5 GeV,z/case) * 6 yA • Az/cos© ). 

If this probability exceeded 0 . 1 , then the step Az was re- 
duced by 2, and the procedure repeated. The cross-section 
of the photonuclear interaction with air atoms 6^ was 
taken constant and equal'vWhen such a photonuclear inter- 
action occured, then the energy of the photon giving rise 
to the nuclear cascade was sampled within the photon spec- 
trum at the depth z, calculated also by meaps of Approxi- 
mation A formulae. The number of 5 GeV muons and the ener- 
gy of the hadron component in that nuclear cascade, rea- 
ching the observation level 700 g/cm 2 , were calculated on 
the base of simulation results, obtained for the "tails” of 
nuclear cascades with the help of the routine for small BAS 
simulation /5/. The result obtained for all partial nuclear 
subcascades in one and the same electromagnetic shower were 
summarized, For 10 13 , 10 14 and 10 15 eV primary energy 1000 
showers were simulated, for 10^-100. 

^-distributions for simulated showers are shown in the 
fig.la. Arrows indicate mean values in these distributions. 
It is necessary to notice that approximation formulae /4/ 
give higher values of Jff e compared with /5 /♦ which differ 
by 1.35 at 10 13 eV and by 1.07 at 10 16 eV. In the fig.la 
results are given, obtained by formulae /4/. It is interes- 
ting to note the change of the distribution shape, when B® 
increases. It is due to the shower maximum, approaching 
the observation level at E®*»10^ ^eV and to the different 
impact of deep penetrations of primary gamma-quantum on 

N , when the shower in the mean is either after the maximum 
6 
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or in the maximum of its 
development* 

The simulation showed 
that the number of photo- 
nuclear interactions with 
energies above 5 GeV is 
proportional to the primary 
energy E* with a good ac- 
Fig. 1 . curacy. The proportionality 

coefficient is o*358 TeV” 1 at 10 13 eV and 0*247 TeV" 1 at 
10 eV,i*e« it is equal to about 0*3 interactions per 
each TeV of the primary energy* The same concerns the frac- 
tion of the total energy, carried out by photons, initia- 
ting nuclear cascades in the shower* It is 0*015 at 10^eV 

and 0*027 at 10 1 eV,i.e* about 2% of the primary energy 
in the average. 



In figures 1b and 1c distributions are shown for the 
number of 5 GeV muons and the energy of the hadron compo- 
nent reaching the observation level^ 700 g/cm 2 . Arrows in- 
dicate again mean values of and B h in the distributions. 
It is seen that the muon number* fluctuates less than the 
energy of the hadron component. In the figure 2 distribu- 


tions of the muon number and of the hadron energy in the 
electromagnetic shower with E®«10 15 eV are shown in com- 
parison with similar distributions for ordinary nuclear 
EAS with the same N & . The latter distributions were obtain- 
ed in the paper /6/. It is seen that and E h distributi- 



ons for electromagnetic and nuclear BAS 
are principally ■"'ather well separated 
from each other* Ratios of their mean 
values, shown in the fig. 3, vary from 
2*0% at 10^ eV up to 4*5% at 10^ eV 
for Njj. and from 1*4# at lO 1 ^ eV up to 
5*0% at 10 eV for E^,i*e* are about 
(2-5)% for both alues* However since 
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the absolute number of muons with the energy above 5 GeV 
In electromagnetic showers is relatively small (~300) 
for E®*10 * eV and they are widely spread over the area, 
then it is usually difficult to identify electromagnetic 
showers by muons. It is possible to make this only at 

arrays with good resolution in the 
region of muon-poor showers, that is 
by means of large area of muon detec- 
tors about several hundreds of m 2 . The 
identification of electromagnetic shower 
by the energy of the hadron component 
is more reliable, since this energy is 
considerably more concentrated around 
the shower core* The calorimeter area 
about several tens of m 2 is enough for 
such identification* 



H. 

Pig. 3 


The simulation showed also that N 


and N e values, as 


well as E k and N 0 don't correlate with each other at the 
fixed value of E«* At the same time there is a strong posi- 
tive correlation between and E h at the fixed E*. 
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